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pH ■ Abstract 

^' 

OO . The Color-Octet Model has been used successfully to analyze many prob- 

lems in heavy quarkonium production. We examine some of the conceptual 
and practical problems of the soft gluon emission process in the Color-Octet 
Model. We use a potential model to describe the initial and final states in 
Cn| ' the soft gluon emission process, as the emission occurs at a late stage after 

^^ ■ the production of the heavy quark pair. It is found in this model that the 

soft gluon Ml transition, ^Sq -^ '^Si , dominates over the El transition. 



^^ ^Pj — > ^5"} , for J/ip and ij)' production. Such a dominance may help re- 

solve the questions of isotropic polarization and color-octet matrix element 
universality in the Color-Octet Model. 

D : PACS number(s): 13.85.Ni, 13.88.-Fe 
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I. INTRODUCTION 

Bodwin, Braaten, and Lepage |]l| have recently developed the Color-Octet Model based 
on nonrelativistic quantum chronio dynamics for very massive quarks which allows a system- 
atic calculation of inclusive heavy quarkonium production cross sections. It is assumed that 
in addition to the production of bound states by the color-singlet mechanism 0], bound 
states are produced by the color-octet mechanism whereby a QQ pair in a color-octet state 
is first formed either by gluon fragmentation or by direct parton reactions, and the color of 
the pair is neutralized by emitting a soft gluon of low energy and momentum. The cross 
section for bound state ip production is represented by 



^(^) = E E nab 



2S+1 TC 



ab c,J,L,S 



(0|O(25+l^c^^^)|0) 



(1) 



where ab = gg,qq,.. are the parton combinations leading to the production, and c = 1,8 
are the color states of the ip precursor. The quantity F{ab —>■ ^'^"'"^Lj ) is the short- distance 
cross section for the production of a QQ pair with quantum numbers '^^^^L'j calculated 
from PQCD, based on the Feynman diagram for ab -^ {cc)cjls + X- The matrix element 
(0|O(^'^"*"^Lj —>■ ^)|0) is the long-distance probability per unit volume for ^■^"'"^Lj to produce 
the final bound state tp and d is the scaling dimension of the operator [|T|. For example, 
for production from the color-singlet state "^^^^Lj to the bound state '4^jlS) which has the 
radial wave function RjLs{f)i the matrix element is 



PI 



^JLst2s+lJj^) 



i^jLsm 



ex 



d'^RjLsir) 



(2) 



For production from the color-octet state '^^'^'^L j, to the color-singlet bound state ipjLSi 
the color-octet matrix element (Og ^"^(^"^ "^^L'j/)) specifies the probability per unit volume 
for the QQ color-octet state ^'^'+^L j, to emit a soft gluon in the transition to the bound 
color-singlet state ijJjls- As the soft gluon emission takes place on a nonperturbative QCD 
time scale and involves nonperturbative QCD, the color-octet matrix elements have been 
treated as phenomenological parameters [p|-[TT|]. Matrix elements have been extracted to 
yield good agreement with the CDF data for high pt heavy quarkonium production in pp 
collisions at y^ =1.8 TeV l^,|l3 and with the fixed-target data for pN collisions at energies 
up to yi = 40 GeV §. 

For the Color-Octet Model to be a valid description, one should be able to describe the 
process of soft gluon emission in physical terms. We know the nature of the final state of 
observed quarkonium, but what is the nature of the initial color-octet state which emits a 
soft gluon? What is the nature of the soft gluon? What is the process of soft gluon emission? 
Which one of the gluon emission processes is more important in J /if) and ip' production? 
Besides these conceptual questions, there are also the following practical questions which 
arise in the application of the Color-Octet Model. 

In the Color-Octet Model for fixed target energies, the dominant production comes from 



lc(8) 



,(8) 



the fusion of two gluons forming a cc color-octet pair in Sq and P) states. The soft-gluon 
radiative transition ^Sq -^ ^SJ leads to an isotropic angular distribution of decay muons in 
the quarkonium rest frame; the ^PJ — > ^SJ transition in contrast preferentially populates 



J^ = ±1 substates with a large transverse polarization, which leads to anisotropic angular 
distributions in the decay of the charmonium to muons. The muon angular distribution can 
be inferred from the magnitudes of the matrix elements (Ogl^So)) and {O^ {^ Pq)) / m^. The 
J/ip and tp' production cross sections at fixed-target energies give the matrix elements @ 



{Oi'^^i'S,)) + ^(03^/^(3po)) = 3.0 X 10-2 GeV^ (3) 

mi 



for J /ip production and 



{Ot\^S,)) + -L(of (3p^)) = 0.5 X 10-2 GeV^ (4) 



for ip' production. On the other hand, the velocity counting rule of ^ gives an {0^{}Sq)) 
of the same order as {0^{^Po))/ml. If we set {0^(}So)) = {0^{^Po))/ml, as suggested by 
the velocity counting rule of [|^ , the expected angular distribution of muons from J/ip decay 
will not be isotropic ^. The experimental data give an isotropic distribution for J/ip and 
ip' production in vr-W collisions at 252 GeV MM and 125 GeV pl. 



For the Color-Octet Model to be a valid description, the color-octet matrix elements 
{OqC^^'^^L'j')) should be independent of the processes which produce '^^'^^L j, . However, 
there are unresolved questions concerning this universality of the color-octet matrix elements. 
High pt CDF measurements at ^/s =1.8 TeV give P 

(Oi^^'So)) + -^{Oi^^ePo)) = 6.6 X 10-2 GeV^ (5) 

m^ 

for J/ip production and 

iOfCSo)) + -^{OfCPo)) = 1.8 X 10-2 GeV (6) 

m^ 

for ip' production. On the other hand, Jf-p and -p' production cross sections at fixed-target 
energies give the matrix elements ^ listed in Eqs. (^ and (Q). To check the universality of 
matrix elements, we need a relation between (O^ipSo)) and {0^(^Po)) for H = {J/'p,%p'}. 
If we again use the velocity counting rule of [|I| to set {O^ipSo)) = {0^{^Po))/ml, the 
fixed-target matrix elements are a factor of 4(7) smaller than the CDF matrix elements for 
J/'p{%p'), as pointed out by Beneke et al. [H. 

We would like to formulate the Color-Octet Model in a form which will allow us to 
answer these conceptual and practical questions. We know that the final observed state 
can be described nonperturbatively in terms of a potential model. Since the emission of 
the soft gluon takes place at a long-time scale after the production of the QQ pair, it is 
reasonable to describe the initial state which emits a soft gluon also in terms of a potential 
model. The soft-gluon emission matrix element can then be evaluated with wave functions 
of the initial and final states in this potential model and the density of color-octet states. 
From these formulations, we find that when a color-octet QQ pair in the ^5*0 or ^PJ 
state emits a very soft gluon to make a transition to the bound H = '^S'} state, the 
probability for the Ml radiative transition ^Sq — > '^S'j is much greater than that for the 
El transition, ^Pj — + ^S*} , and the velocity counting rule breaks down. We suggest that 



such a dominance of (Cg ^(^5'o)) over (Cg ^{^Pj))/ml for soft gluon radiation may help 
resolve the above questions of quarkonium polarization and matrix element universality. 

This paper is organized as follows. In Section II, we formulate the model of quarkonium 
production in terms of the short-distance Feynman amplitude and wave functions determined 
from potential models. An explicit expression is obtained to relate the color-octet matrix 
element in the Color-Octet Model to quantities which can be evaluated in the potential 
model. In Section III, we examine the potential between Q and Q interacting at large 
distances. When one allows for the effect of the spontaneous production of light quarks to 
break up the string joining Q and Q, the linear potential is modified to become a screened 
potential. We use a screened potential to discuss bound states in Section III and continuum 
resonance states in Section IV. Section IV also outlines how we obtain continuum wave 
functions which are needed in the calculation of the cross section. Section V gives the 
multipole transition matrix elements. As an illustration, a simple Ml to El ratio is obtained 
for the production of J/ip using approximate wave functions. Numerical results with more 
realistic potentials are obtained in Section VI to give the contributions from Ml to El 
transitions to J/tp and ip' production. It is found that the soft gluon Ml radiative transition 
^5*0 — >■ ^5*} dominates over the El transition ^Pj -^ ^Si for J/ip and ^' production. 
Section VII summarizes and concludes the present discussions. 

II. HEAVY QUARKONIUM PRODUCTION FROM PARTON COLLISIONS 

To obtain the cross section for the production of quarkonium bound states in a hadron- 
hadron collision, it suffices to focus on the production cross section in parton collisions, as 
the former can be obtained from the latter by folding the parton distribution functions of 
the colliding hadrons. We consider the collision of the parton a with the parton b which 
form the initial QQ pair. The heavy quarkonium production amplitude can be obtained by 
projecting out the state vector of the QQ pair onto the quarkonium bound state. We shall 
show in detail how this is carried out to yield the quarkonium production cross section. 

The production of a heavy quark pair is a fast process. We can follow the time evolution 
of the state vector of the QQ pair. The state vector resulting from the collision of a and b 
at initial production time ti is 



\^a,{u)) = \^l'^m + \^Tm + ... (7) 



where 



Here, Pi and Pf are the initial and final 4-momentum. The matrix elements Ai can be 
obtained in perturbative QCD using Feynman diagrams and renomalization procedures. 
For a sufficiently large value of t, the initial state in Eq. (0) will evolve to become 



\^ab{t)) = E fd^^^- (27r)^5(^)(P, - Pf)M{ab -. tjjjLs x)\tljjLs x), (10) 

where x denotes the number of hard gluons, and d^^x is the corresponding phase space 
volume element. For example, for ab ^ ip, d^^ is 

and for ab ^ tp g with the emission of a hard gluon, d^^g is 

"^ ^'~ (27r)32E^, (27r)32E/ ^ > 

The cross section for the production of the quarkonium state ipjLS by the collision of partons 
a and b is then 

da{ab -. ijjLs x) = jj{2rry6(^\n - Pj) E \M{ab -^ tfjjLs x)\^d^^x (13) 

X 

where / = JiPa ■ PbY — Tnlml. It is necessary to obtain A4{ab -^ ipjis x) from Ai{ab — >• 
QQ x) in order to calculate the production cross section. We shall show how this relation 
can be obtained in the present model description. 

One can first separate out the center-of-mass motion and the relative motion of the QQ 
pair in terms of their total momentum P and their relative momentum q. The heavy quark 
momentum is then Q = P/2 + q, the antiquark momentum is Q = -P/2 — q, and the state 
of the center-of-mass motion can be represented by the plane wave state \QQ) which is 
equivalent to \qP). They are related by the ratio of their normalization constants. 

We look first at the situation when Q and Q are not interacting to introduce their 
interaction through a mutual potential V later on. Part of the phase space elements, d'^Qd'^Q, 
in Eqs. (|^) and (^ can be transformed into 

^'^-5(g^ - mD^AQ' - ml) = ^MP' - ml)^J (^) (14) 



(27r)3^^ ^^(27r)3^^ «' (27r)3^ ^' (27r) 

where rrip is the invariant mass of the QQ system. To study the dynamics of the relative 
momentum, it is best to refer to the center-of-mass system where P = (mp,0), q = (0,q), 
and the relative momentum q satisfies the following equation 



where 



and m^ is the reduced mass. 



C - g - < = 0, (15) 



— ^ ^, 16 

2mp ' ^ ' 



m^ = niQ/mp. (17) 



Next, we consider the system with a mutual time-hke vector interaction V. The equation 
of motion can be obtained from Eq. ([T5|) by the canonical method of replacing e^ with 
^w ~ V{r). The Klein-Gordon equation for the two-particle system under a mutual time-like 
vector interaction V{r) is 



[e. - V{r)Y -q'- mij<if{r) = 0. (18) 

Similar two-body equations of constraint dynamics for more complicated cases with spin 
and very general types of interactions can be found in [jl^ . Because of the large mass of the 



heavy quark, it is convenient to use nonrelativistic QQ wave functions and state vectors. 
The Klein-Gordon equation becomes the Schrodinger equation, 

where e = e^— m^^ is the non-relativistic measure of energy, and m^ ~ Mq/2. The eigenvalue 
e of the bound state is related to the masses by mp = e+2mQ. The bound state wave function 
can be written in the form 

^jLs{r) = RjLsir)yjLsir). (20) 

The wave function in momentum space is then 

i^JLs{q)= fdre-''l-''^jLs{r). (21) 



Heavy-quarkonium bound states have been obtained previously with many different forms 
of the potentials [|l^ . 



In nonrelativistic kinematics, P = {Pq, P), and q = (go, q) ~ (0, q). One can perform a 
decomposition of the relative wave function in terms of color and angular momentum states. 
The amplitudes in Eq. (|^) can be decomposed as 

M{ab ^ QQx)\QQ) = Y. M'j^sil ^)\<1 P), (22) 

cJLS 

where M'jisi.Q ^) is 

M'j.siq x) = A'j,si\q\ x)y^jLs{Q) (23) 

with yjLsi^) ^^^ angular momentum and color wave function. For simplicity of notation, the 
azimuthal angular momentum component M and the color component of the color multiplet 
in yji^s ^ill ^ot be written out explicitly. The bound state ipjis can be described as 

li^JLS P) = I ^3^JLsiq)\q P). (24) 

In the lowest order of approximation without soft gluon emission, the probability amplitude 
for the direct production of 'ipjLS is obtained by projecting iIjjls onto (^abiti)- We use the 
normalization for the center-of-mass motion 

6 



(P'|P) = {2nf2Ep5{P' - P) (25) 

and the usual convention for nonrelativistic relative motion 

{q'\q) = {27cr5{q'-q). (26) 

To carry out the projection of ^abiti) onto ipjLS: we use Eq. (|2^ ) to write Eqs. (H) and (^ 
as 

The scalar product of \ipjLS Px) with the above state vector, {ipjis PA^ahi^i))i '^^'^ t)e 
evaluated by transforming to the variables P and g with Eq. (|T^. Using Eq. (P^, we obtain 

(V^,i5 Px|<|.««^(t,)) = {27:Y6^'\P^, ~^^^^sl -§^^'i'*J^s{q)Mf,s{Q ^)- (29) 

Therefore, from the above result and the definition of M.{ab —>■ ipjis x) as given by Eq. 
(0), we get 

Miab ^ ^jLs x) = —{^jLsiq)\M^Jlsiq x)), (30) 

where 

{^JLs{q)\M%{q x)) = J -0-^rjLsiq)M%{q x), (31) 

which involves only color-singlet components of \^ab{ti))- Thus, if the bound state wave 
function is known, the above overlap will give the matrix element and the color-singlet 
contribution to the quarkonium production cross section in Eq. (|l^). 

In many approximate calculations one expands Ajlg{\q\ x) in powers of the velocity 
V = \q\/MQ and retains only the leading term. 






where 

Mk r ^L ,^^ 1 

(33) 






^JLSx 



L\ 



"^^ A^IUWI ^) 



d\q 



\q\^o 
Then the above projection in Eqs. ( pOD and ( pi]) gives the L-th derivative of the wave function 



at the origin [|T9|-|2T 



{^JLs{q)\M%{q x))/M'q^ = Af,sA^jLs{q) \ l^l'yjLS 



(1) .. (2L + 1)!! 



Mc 



dr^ 



(34) 



Jr^O 



The cross section can then be separated into a short-distance part involving Ajj^g^ and the 
long-distance nonperturbative part involving the L-th derivative of the bound state wave 
function at the origin. We can write the above results in the notation of the Color-Octet 
Model of Eq. (g): 

.(,) . 2 2 na» ^ -L5 ^m^l^I^^lM, ,35, 

ab c,J,L,S ■'^-'Q 



where 



na, ^ -«L« ) ^ / 1(2,)'..'.(P, - P,) A E .XT:'!)!! "**- P*" 



and following Bodwin et a/. |jl| 



(010(^^+^4'^ ^ V; )|0) _ 2Ar, (2L + 1)!! 
M2^ ~ 47r L\ 



Mkdr^ 



(37) 



r^O 



In addition to these color-singlet contributions, the color-octet model further postulates 
a much slower process with the emission of a soft gluon Qs-, different from the emission of 
the hard gluon g in perturbative QCD already included in the second amplitude |$^j, ) of 
Eq. (0). One can describe the emission of a soft gluon by studying further evolution of the 
state \^abiii)) under the influence of the color field in the form 

\^ab{t)) = U{t,U)\^ab{ti))--{l-ig j^ Jj-Ad^xd?! l^abiU)), (38) 

where we shall focus attention on soft gluons only for the gluon field A. 

The probability amplitude M{ab — > ipjLs Qs) for the production of the bound state 
ip{JLS) accompanied by a soft gluon Qs is 

{2T^Y5{Pab -P^- Pg)M{ab -^ ^jLs gs) = {[ipjLs; P^] 9Ab{t)) 

= {[^pjLs;Pq] 9s\ {-ig) j J ■ A d^x dt l^abiU)) (39) 

where we have used the notation Qs with the subscript s to denote a soft gluon emitted by a 
color-octet state. Using Eqs. (H) and ([l^), the righthand side of the above equation becomes 



RHS of (39) ^ E lj^,HP' - "Djl^* (^) (2^)'*'"(P. - P) 



J'L'S'x' 



(27r)3^ ^^(2 



^{{.il^jLsP^)gs\{-ig)JlJ''3 ■ A d^x dtM%s'iQ x)\qP)- (40) 



One envisages that in order to emit a soft gluon, the state $ab must have evolved for 
a long time and has characteristics of nonperturbative states. It is therefore reasonable to 
postulate that the soft gluon is emitted by an intermediate nonperturbative state Wj'l's' t'^st 
described nonperturbatively in terms of Q and Q interacting in a potential. It is appropriate 



(8) 



to expand the intermediate states in terms of V'j'l's' with an invariant mass mj'L's' using 






(41) 



Here dn/dmj'Lis' is the density of these color-octet wjiys' states with a specific polarization 
and color-octet component: 



dn 



m, 



3/2 



dmjiL's' (2vr)2\ 



mj>L'S' - 2mQ - 



L'{L' + 1) 
mqE? 



(42) 



where R is the average radius of the Q-Q separation of the color-octet state which emits the 
soft gluon. It can be taken approximately as the Q-Q separation of the final color-singlet 
quarkonium state. For the emission of a soft gluon through the intermediate color-octet 
state V'j'L'S'; th^ righthand side of Eq. (|39|) is 



RHSof(39)= E j^,S{P 



— nij, 



d^q f P ■ q 



J'L'S' 



J'L'S'x • 



(27r) 



tS 



{2T:)%Pab - P) 



dn 



X {{'il^jLS P^)gs\{-ig) / J -Ad^x dt\ij%s' P)dmj,L's' , 



,(8) 



,(8) 



{r;/L's'\M'j>,,s'i<i^)\<i)- (43) 



The integral over q can be carried out to give 
d'^q f P ■ g' 



(2vr) 



:S 



{^%s'\M%s'i<l^)\<l) 



d^q 



',(8) 



i'%s'{<i)M%s'{<i^) 



,(8) 



mj'L'S' 



{rA'S'^M^P^.s^Mx)). (44) 



Thus, the amplitude involves a projection of the color-octet Feynman amplitude onto a 
color-octet nonperturbative amplitude, {i^jiLis'i^)\-^j'L'S'i^^))i similar to the projection 
of the color-singlet component onto a color-singlet amplitude as given in Eq. (|30|) . One can 
thus similarly relate this projection to the L-th radial derivative of the spatial wave function, 
as was shown in Eq. (^4]). One obtains 



.,(8) 



,(8) 



(V^:7'i'5'(?)l-^>'W(? ^)) _ ^(8) , .,L' (2^' + 1)!! 



M^ 



MK^-nL'\ 



d^ Rj>L'S'ir) 
dr^' 



(45) 



J r-^O 



where 



•^J'L'S'x 



L'\ 



l(8) 



d\q\ 



tA'L'S'{\q\ ^) 



(46) 



\q\^o 



To evaluate the long-distance soft gluon radiation matrix element, we have 



Ait 



fc V '^^9 

where ujg is the soft gluon energy, e^ is the gluon polarization vector, Ck and c[ are the 
creation and annihilation operator for the soft gluon. Therefore, we have 

{{^jLsP^) gs{Pg)\ {-ig) f j-A d\ dt \^%s'i(i) P) 

J ti 

= Vj,L's'^JLs{2nfSiP -P^~ Pg) f e-^^^^-^f-^^^'dt 

J ti 

where 



VrL'S'^jLs = -^^y ^< Jfj'L's'-^jLsir)e''^-''dr, (48) 

the time-like component oi j%{r) = Pfi{r) is 

Pfi{r) = / dri dr2ip}{ri, rs) Y. ^i'^ - rn)4^i{ri, ra) (49) 

•' n=l 

where ri,r2 refer to the coordinates of the heavy quark and antiquark, and ipi and ipj are 
the initial and final wave functions. The space-like component is 

Jdr) = -^{rf{r)VUr) - ^r(r)V^;(r)) + |^V x (^;(r)^^,(r)), (50) 

where /ig ~ 2 is the color-magnetic moment of the quark in units of the color Bohr magneton. 
Therefore, we have 



M{ah -^ VA^S'^ -^ i^JLS Qsx) = {ijyL'S'iQ^)\^j'L's'iflWj'L'S'^JLS 2 27r-- (51) 



where Ejilis' = y "^ja's' + (« + b^- From Eq. ([T3|) 

da{ab^^jLs9sx) = ^ ^ \Miab ^ ^%s'^ ^ ^JLs 9sx)\^i27T)''5^^\P, - Pf)d^^g^, (52) 

^^ xJ'L'S' 

where the phase space volume element d^^g^x includes both the soft and the hard gluons. 
For the case when only a soft gluon is emitted, d^^g^ is given by 

^*^^^ - (2vr)32i?, (2^)32E,; ^^^^ 

We have then 

da{ah ^ i,j,s 9s) = ^. E (^— )'l(^i'U'(g)l-^S'U(g^))r 

4^ xJ'L'S' ^^J'L'S'y 



y<\Vj'L's'^jLs\'-r^ — {27T—^^^—\ ^dn^. (54) 
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where \p^\ is the magnitude of the momentum of iI^jls in the rest frame of ipjiL's'- O^^ ^^^ 
expand {'4'j'L's'i^)\-^j'L'S'i9 ^)) "^^^ according to Eq. ( ^3D and obtain 



d(T{ab -^ i/jjLs Qs 



E 



^^ xJ'L'S' ^^J'L'S 



x\Vj. 



L'S'^JLS] 




r^O 



8m j,^, 5, \dmjiL 



(55) 



We can then compare the above results with the expressions in the Color-Octet Model, 
which we write in the form 



da{ a b-^ ipjLs 9s) 



1 E na„ ^ .rL'S')i2.fm,, + p. - p,) ^^^*^^ ^^''™ 



where 



AT ^^ 
^-' J'L'S' 



fiab - J'L'S') = 1 E (-^^ ' 



/ (2L- + l)!! -(8) 



(56) 



(57) 



where Cj'l's' = (number of polarization) x (A^^ — 1) for the color-octet state {J'L'S'}. We 
observe the following differences. The Color-Octet Model result of Eq. (|56D assumes that 
the soft photon does not carry energy, hence the cross section is a delta function of the total 
parton invariant mass, centered at the invariant mass of ipjis- The soft gluon makes no 
contribution to the phase space volume element and the matrix element is parameterized 
as unknown constants. On the other hand, when one takes into account the details of the 
emission process to obtain the results in Eq. (|55|) , the cross section is a continuous function of 
the total parton invariant mass mjiiis'- The production of ipjis occurs within the invariant 
mass interval from the minimum ?7ij//,/5/(min), at which the color-octet state begins to exist, 
to the maximum value mj/L/5/(max)=mji;,5'-|-ci;g(cutoff), which corresponds to the maximum 
gluon energy of the soft gluon process. One envisages that gluons with energies ujg higher 
than cUg (cutoff) will be emitted at a shorter time scale and will be described by perturbative 
QCD and not by the soft gluon emission process through an intermediary non-perturbative 
state. As the time to travel the length of R/2 of a J /if) corresponds to a energy scale of 2h/R, 
this tUg (cutoff) energy should be smaller for a final quarkonium state with a greater radius 
(such as ip'). Furthermore, another difference is that the cross section of Eqs. ( ^2]) and (|55D 
in the present description is given in terms of transition matrix elements of nonperturbative 
wave functions, and the phase-space volume element includes the effect of the additional 
soft gluon. 

In an approximate comparison by integrating Eqs. (|55|) and (|56|), we get approximate 
equivalence 



{0^^\J'L'S' -^ JLS)) 



M'q'' 



E 



JLS 



"Cj'l's' (™ax) 



dm 



I (min) 



J'L'S' 



ICj>L's'{2L' + r 



A-kL'I 



d^R 



J'L'S' 



(r) 



M^'dr^' 



x\Vj'L'S'- 



dn 



>JLS\ 



^T^rnj,^,g, \draji 



US' 



iPv 



(if2„ 



(5^ 
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We shall consider the soft gluon to have large wave lengths which are greater than the 
dimension of the radiating system. In this long wave length approximation the matrix 
element Vfi has been worked out in detail (See Eqs. (46.6) and (46.7) of Ref. |^). The 
results for gluon radiation of multipolarity LM are 

\Vf^\ -4 ^ [(2L + l)!!p'^^^ ' ■ (^^^ 

For electric transition of multipolarity LM, the multipole transition moment is 

^Sf = gJdrpf,iryYLM{r)- (60) 

For magnetic transition of multipolarity LM, the multipole transition moment is 



r(m) 



l—Jdrrx j,,{r) ■ V {r'-Yu,(i-)) ■ (61) 



M,„ ^^ 



III. POTENTIAL MODEL FOR Q AND Q INTERACTING AT LARGE 

DISTANCES 

We wish to describe the nonperturbative interaction between Q and Q in terms of a 
potential model. Such an interaction is known for Q and Q in the color-singlet state. 
The interaction consists of an attractive short-range color-Coulomb interaction and a long- 
range linear potential. For Q and Q in a color-octet state, the short-range color-Coulomb 
interaction is repulsive. From lattice calculations for gluonic excitation in the presence of a 
static quark-antiquark pair, the spectrum appears to rise linearly with the relative separation 
between Q and Q, indicating an effective linear potential for color-octet states |l^. We can 
write the Q-Q potential as 

V,ir) = C,y + a,r (c=l,8) (62) 

where Ci = —4/3, Cg = 1/6, a^ is the strong interaction coupling constant, ai ~ 1 GeV/fm, 
but as is not known. 

The confining potential given in Eq. (|62D is a good description if the spontaneous pro- 
duction of light quark pairs is not considered. However, when a quark and an antiquark 
pulls far apart, the string joining Q and Q breaks. Evidence of string breaking comes exper- 
imentally in the production of a large number of pions in e^e~ annihilation at high energies. 
This phenomenon can be explained as follows. The annihilation of e"*" with e~ produces a 
qq pair. As the q and the q pull apart, the string between the quark and the antiquark 
breaks up into sections of strings which manifest themselves as pions p3|j2^ . Theoretically, 
the pair production will occur spontaneously through the Schwinger mechanism when Q is 
separated from Q at a. distance r such that ar is greater than the mass of a qq pair |^ . 



When the string breaks, the interaction between Q and Q will be screened and will become 
a weak color-van-del- Waal-type interaction, instead of the strong linear interaction. The 
asymptotic behavior of the Q-Q potential is therefore drastically changed when the sponta- 
neous production of light quark pairs is allowed. Instead of the confining behavior as usually 
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assumed, it is more appropriate to describe the motion of the Q and Q effectively as having 
an asymptotically free motion at r — ;> oo, on account of the breakup of the string due to the 
spontaneous production of light quark pairs, leading to the production of a D and a D. 

The influence of the dynamical light quarks on the Q-Q potential has been investigated 
theoretically by Laermann and collaborators [0,0 using lattice gauge theory. It is found 



that including the dynamical quarks in the lattice gauge theory leads to a potential which 
deviates from the linear potential. The modified interaction can be described in the form 
of a screened color- Yukawa potential. Because of this screening due to dynamical quarks, a 
more realistic Q-Q potential becomes EO, 



V;(r) = C,^^^ + a/-^ ^, (63) 

r jicT IJ.C 

where we have used the reference that the potential vanishes at r — i> oo. 

Under the potential V^(r) given above, nonperturbative QQ states ipjLsi''^) can be ob- 
tained by solving the Schrodinger equation 



V' + V,ir)\rjLsir) = erjLsir), (64) 



where e = rrijis ~ 2Md. 

In our approximate treatment, we consider a single heavy object which can travel at 
all distances, with the provision that in going from small distances to large distances, the 
nature of this object changes from being a charm quark at small distances to becoming an 
open charm D-meson at large distances. Thus, the mass of this heavy object should also 
change from that of the charm quark mass to a D-meson mass. We describe such a variation 
of the mass by the effective mass 

where we shall use r^ = 0.8 fm and am = 0.2 fm to describe the change in nature of this 
heavy object at a distance of r^ ~ 0.8 fm. 

We use such the potential (|63|) and the effective mass (|65|) to study the charmonium 



bound states and the resonances in the continuum. The latter are characterized by an 
asymptotic phase shift of (2n + l)vr/2 where n is an integer. A reasonable description of the 
bound states and resonances can be obtained (Table I) by using the parameters of a^ = 0.3, 
(Ji = 1.7 GeV/fm, /i = 0.28 GeV, and Mc=1.6 GeV. A potential similar to (|63|) has been 



used successfully to study the energy spectrum of charmonium states |]2B|. As in this earlier 
work, we find that a large value of the parameter cxi is needed when the linear interaction is 
modified to take the form of the screened potential of Eq. ( p5D . We shall use this potential 
to generate the J/ip and ip' wave functions for subsequent studies. 



Table I. Charmonium states obtained with the screened potential Eq. 
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Bound State 


Resonance 


Calculated 


Experimental 


Calculated 


Experimental 






Energy (GeV) 


Energy (GeV) 


Tee+ (keV) 


Tee+ (keV) 


IS 




3.145 


3.07 


8.11 


5.26±0.37 


2S 




3.560 


3.66 


2.57 


2.14±0.21 




3S 


3.947 


4.04 






IP 


2P 


3.476 
3.971 


3.52 








ID 


3.801 


3.77 







The behavior of the color-octet {QQ)8 states can be studied by looking at the states in 
the potential (|63|). The color-octet potential is repulsive at short distances. For sufficiently 
large values of as, there is a potential pocket at intermediate distances. Depending on the 
mass of the quark and the string tension ag, the pocket of the color-octet potential may hold 
a bound state. The greater the rest mass of the quark, the greater the chance of a color-octet 
bound state. In this respect, there may be a greater chance of finding a color-octet bound 
state in a bb system compared to a cc system. 

If the linear potential parameter as for color-octet states is greater than 0.31 GeV/fm, 
then the potential pocket will hold a bound [cc\s state. (The state is bound in the sense 
that its mass is less than 2M£), but it is unstable against the emission of a gluon.) There 
is, however, little experimental information on the color-octet states of either a cc or a 
bb system. Compared to color-octet states in the continuum, a bound color-octet state is 
distinguished by its relatively longer lifetime, as the only decay mode is the emission of soft 
gluons, and the emission rate is lower for the bound state because the energy difference 
between the initial and the final state is less for a bound color-octet state. A long-lived 
color-octet quarkonium will travel a large distance before its color is neutralized. 

It is of interest to examine the signature of a long-lived bound color octet state. Its only 
decay channel is the emission of a gluon to come down to color-singlet QQ bound states. 
The emitted gluon will be converted to a light qq pair. If the initial color-octet cc pair has 
a high kinetic energy relative to its complementary color-octet partner (which is either a 
[qq]s or a [q{qq)]8), then the qq from the emitted gluon will form two chains of light mesons 
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(mostly pions) with the color-octet partner as shown in Fig. 1. 

vVvvvV/ 



71 or A? 71 7t 71 71 



5 or diquark 




[q(q or diquark) ]g 



[QQh 



Fig. 1. Schematic picture to show how a bound color-octet [QQJs and its 
color-octet partner [ggjg (or [q{qq)]8) evolve into color-singlet objects. 

In the more general case when the color-octet is accompanied by two color-octet partners 
(one from the projectile hadron and one from the target hadron), the emitted gluon from 
the [QQls color-octet needs to branch out into two gluons, which will be converted into two 
pairs of light [qq]s pairs to form four chains of produced light mesons. 

To search for long-lived color-octet states, it is useful to look at the kinematic regions 
where the produced heavy quarkonium is kinematically separated from the other debris of 
the collision, as in the case of J/ip or T at large values of xp. The signature of a long-lived 
color-octet quarkonium state may show up as a quarkonium at one end accompanied by 
hadrons which arise from the breaking of strings joining this color-octet and its comple- 
mentary color object at the other end. There will be chains of pions in the rapidity gap 
between the produced quarkonium state and its complementary partner. Thus, the search 
and the measurement of the chains of hadrons with a leading quarkonium state may re- 
veal the existence and the characteristics of the quarkonium as a color-octet state. In the 
situation that the pocket is not deep enough to hold a bound color-octet state, then the 
intermediate color-octet state ^j/i/^/, through which the soft gluon is emitted, can only be 
in the continuum. 



IV. CONTINUUM STATES 

For an intermediate color-octet state in the continuum, the cross section for the emission 
of soft gluons depends on its wave function in two ways. First, it is needed to give the overlap 
between the amplitude from perturbative QCD as given by Eq. (|^). It is also needed to 
obtain the transition matrix element for soft gluon emission (given in the next section). 
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To obtain the continuum wave function, we follow the phase-angle method discussed in 



detail by Calogero [29|. We write the wave function as 



i'jLsir) = RjLs{r)yjLs{r) = "^^^I^yj^s^r) (66) 

and represent the wave function in terms of the amplitude aL{r) and the phase shift S^^r) 



ujLs{r) = ^f^Mkr) sm{6L{kr) + h{r))J^, (67) 

with the boundary condition that 5l{t — » 0) = 0. The functions D{kr) and S{kr) are known 
functions [p9| : 



Doix) = 1, Di{x) = (1 + l/x^Y^^, (68) 

and 

Sq{x) = x; Si{x) = X — tan~^a;. (69) 

The equation for Sii^r) is [^ 



^Sr.{r) = -^lDl(kr) {sinMkr) + (^^(r)]}' , (70) 

where U{r) = MQ{r)Vc{r) + (M^ — MQ{r))e. After the function ^^(r) is evaluated, the 
amplitude can be obtained from ^^(r) by 

aL{r) = exp 1^ J^ ds U{s) Dl{ks) sm2[6L{ks) + Sl{s)]'^ • (71) 



The numerical integration of Eq. ( [fOD gives the asymptotic phase shift 5{oo) and the 
continuum wave function. The energies at which the asymptotic phase shifts 6{oo) are 
(2xinteger+l)7r/2 are the locations of the resonances as shown in Table I. 

For the wave function and its derivatives at the origin, it is useful to express the effect of 
the potential in terms of the i^-factor, the analogue of the Gamow factor for the Coulomb 
interaction |^. The wave function [Eq. (|67D ] near the vicinity of the origin is given by 



RjLsir) = — = — -, — rjL{kr^ 



kr aiioo) ' ^ 2L + 1 



aL(oo)(2L + l)!! V 2L + 1' 
Therefore, for the factor in Eq. (|55|), we have 

/(2L + 1)!! 1 d'^RjLsir 



(72) 



47rL! M^ dr^ 



«L(r) k^ / L\ 

r^o ~ M^)MSn'2L + m{2L + l) 



(73) 



Compared to the case of no interaction for which a{r) = 1, the probability is modified by a 
factor, the /('-factor, 
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Kr 



aUr -^ 0) 



(74) 



In terms of the K-factor, the square of the L-th derivative of the wave function at the origin 
is 



'(2L + 1)!! 1 d'^RjLsir) 



An L\ Mk dr^ 



r^O 



Kr 



k 



M, 



Q 



2L 



L\ 



(2L + 1)!!(2L + 1)' 



(75) 



An example of the /('-factor is the Gamow factor arising from the Coulomb interaction. 
Our numerical i^-factors can be checked by noting that for the Coulomb interaction with 
V{r) = a/r, the i^-factor can be obtained analytically and shown to be 



Kl{v) 
where t] = a /v. 



{L^ + r]')[{L-lf + r]^]...{l + 7]^) 



2TTri 



[L\{2L + 1)\ 



112 



1 — exp{—2Tiri} 



(76) 



V. MULTIPOLE TRANSITIONS 

We can evaluate the square of the matrix element | Vfi \ ^ for the emission of soft gluons 
based on the analogous multipole expansion for the emission of electromagnetic waves p2| , |3T 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

r (fm) 



Fig. 2. Schematic picture of the color-singlet potential V^^\r), color-octet 
potential V^^\r), the initial color-octet state Wjilis'^ ^^"^ ^^^ final color-singlet 
state ijjjls i^ ^ soft-gluon emission process. 
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We consider the production of a QQ pair by parton collisions at an energy e above 
the DD threshold. The QQ can be in many different angular momentum states. After 
a nonperturbative time scale, the QQ pair will evolve into nonperturbative states ipj/L/s'- 
We consider the continuum color-octet states V'j/iQ 2,/=o s'=o ("'^'^'o ) and ■j/'j; ^.'=1 s'=i i^-^j ) 
which can decay down to the color-singlet ■i/'jii ^.=0 5=1 (^'^'i ) state by the emission of a 
gluon of energy cUg = e + B where B is the binding energy of '^5'} . (For brevity of notation, 
the color state superscript in the wave function ipj^s is often understood.) The square of 
the matrix elements for the emission of a gluon of energy Ug = e + B from the continuum 
color-octet QQ(^Sq ) state to the bound color-singlet H={^Si ,Sz = 0) state through an 
Ml transition is given by [0,0] 



p^Mi^i^W ^ 3^(1)) - ^ ( ^ )ujI\M[,\\ (77) 



.9 
The magnetic dipole transition matrix element M[q is 

^^°= 2Iq\/? E_/^Hq('^o'V.cx..^h(^^;^^), (78) 

where /ig is the color-magnetic dipole moment of Q and /ig = — /iq ~ 2. The ipQQ wave 
function for the color-octet QQ{J' = 0,L' = 0, S' = 0) pair in the continuum can be written 

as 

^I^qq{J'L'S') = ^^^°^ ' ' yj'L's'. (79) 

The wave function for the bound ipjis state can be written as 

i^JLS = ^^^^^^yjLS. (80) 

r 

From Eqs. ([77|)-(pOD, the matrix element M(q is given by 

'uCSo'^^\r)uCs[^\r)dr. (81) 

We shall evaluate the above matrix element with wave functions obtained from the potential 
model. As an illustration of the order of magnitude, we can consider the simple case of 




uCSo^^\r) = sm{kr)V^ (82) 

and the lowest bound ^Si state wave function as 

«(35/'\r) 



2VK^e-^' (83) 



where k = JMqB. Then the square of the matrix element for Ml gluon radiation from 
'S^'^ to 3S« is 



^Ml^l^(8) 



3^(l)^ 



1%^(/^Q - f^q) 



2, ,3 



K 



M'r 



(P + k2)4' 



(84) 



The Ml transition rate is constant as fc ^ because it involves the spatial overlap of two 
S'-state wave functions. 

The square of the matrix element for the emission of a gluon of energy ujg from the 
continuum [QQi^Pj '), J^ = Sz) state to the bound if=(^5'j , S^) state through an El{m = 
0) transition is 



El{3-ni8) _^ 3q(1) 



v'^'rp 



J 



SI' 



U^alQ 



10 1 



(85) 



The electric dipole transition matrix element is 



Qio = W£ T._Jd'r^hQ("Pj'^) '^ *^('^f 



3q(1)n 



ig_ 3_ 
k\4n 



/3p(8) 



u{^Py',r) ru{'Sl'\r)dr^ 



3q(1) 



(86) 



As an illustration, we shall again evaluate this matrix element with a wave function for the 



non- interacting case where the Ugp(s) is 



ueP!f\r) 



'sinA;r \ Air 

— ; coskr \ — . 

kr / V 3 



(87) 



The square of the matrix element for El gluon radiation from ^Pj to ^Si is 



^1/3 0(8) ^ 3q(1) 



^i.l/3p 



J 



SI' 



4 X 256^2 ^ ^5^2 
27 ""^(FT^' 



(88) 



The El transition rate vanishes as k ^ 0. This arises because the transition matrix element 
involves the product of the S-state wave function, which peaks at the origin, with the 
continuum P-state wave function, which vanishes at the origin. The product is zero in the 
long wavelength (fc— >0) limit and increases with increasing kinetic energy. 

For the simple case when there is no interaction in the color-octet state, the ratio of the 
squares of the matrix elements for the production of the lowest tp state are 



Afl/l c(8) . 3cW 



V^'H'So 



SY 



Y El ^3 pip 



■igw'' 



16 eMg' 



(89) 



With hq = 2, we see that this Ml /El ratio is large for soft gluons which carry energies 
slightly greater than the binding energy B. The Ml transition arises from the spin transition 
current. In contrast, the Ml transition arising from the spatial transition currents would 
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lead to an Ml/El ratio varying as f^ = Aq'^/Mq. In the present case, because of the angular 
momentum and spin selection rules, there is only a contribution to the Ml transition from the 
spin transition current, and the usual velocity counting rule based on the spatial transition 



currents breaks down 31 



The ratio of the | V^jp for the production of the ip' state will be more complicated as the 
ip' wave function contains one node and the radial matrix elements will involve cancellation 
of contributions of opposite signs. 



VI. CONTRIBUTIONS FROM Ml TO El TRANSITIONS TO J/^j AND V' 

PRODUCTION 



From Eq. (pB|), we note that the cross section for going through different intermediate 
color-octet states J'L'S' depends mainly on the product 



C 



J'L'S' 



'(2L' + 1)!! 1 
47rL'! Mk' 



d^ R J'L'S' jr) 
dr^' 



Wv 



dn 



J'L'S'^JLSl 



drriji 



US' 



(90) 



We can consider the production of the J/ipi^Si) state via an Ml transition from the in- 
termediate color-octet state {J'L'S'} = ^Sq. There can also be an El transition from the 
intermediate color-octet state {J'L'S'} = ^Pq. To evaluate various quantities in the above 
equation, we need a description of the color-octet state. There is the short-range repulsive 
interaction represented by Cc = —1/6 in Eq. (|63D. If the linear interaction between the Q 
and Q is strong enough, a bound color-octet state may be present in the pocket in the in- 
termediate distance range to feed to the observed color-singlet state. The bound color-octet 
state will most likely be an L' = state rather than L' = 1 state because the latter lies 
at a much higher energy. A bound color-octet L' = state will make the transition to the 
color-singlet state by an Ml transition, and there will be no El transition in this case. 

Due to the absence of experimental evidence for bound color-octet states at present, we 
shall consider color-octet states to exist only in the continuum and shall assume only the 
color- Yukawa interaction with no linear interaction. Using such a color-octet potential, we 
obtain the wave function in the continuum using methods outlined in Section |ITI[ The K- 
factor can be evaluated to give the square of the L'-th derivative of the wave function of the 
continuum color-octet state at the origin [Eq. (|75D]. The wave functions in the initial color- 
octet state and the final color-singlet state can be used to evaluate the matrix elements of 
Ml and El operators. We then obtain the various quantities in Eq. (|90|) . For the purposes 
of displaying the results, we introduce the product Fr defined as 



'R 



C 



J'L'S' 



'(2L' + 1)!! 1 
47r L'\ Mh 



d^'Rj,L's'{r) 
dr^' 



(91) 



There is a limit for the maximum energy of an emitted soft gluon. One envisages that 
gluons with an energy higher than u;^ (cutoff) will be emitted at a shorter time scale and 
will be described by perturbative QCD and not by the soft gluon emission through an 
intermediary non-perturbative state. The cutoff of this soft gluon energy depends on the 
size of the bound state, as the time required for a gluon to travel the length of the radius 
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R/2 of a J /ip corresponds to an energy of 2h/R ~ 0.8 GeV. (Here, R is the separation 
between Q and Q which is about 0.5 fm for J/ip.) Therefore the maximum energy of an 
emitted soft gluon should be of the order of 0.8 GeV. The extracted color-octet matrix 
element depends on this cutoff gluon energy. We find that a value of ujg{m.a.x) = 0.7 GeV 
gives the best description and is approximately consistent with this estimate of 2h/R from 
the size consideration. For this maximum value of soft-gluon energy, the invariant mass of 
the color-octet state is 3.8 GeV. 

In Fig. 3 we show |l^jP, -Fr, and {dn/dmY as a function of the invariant mass mj/irs' 
of the intermediate color-octet state for the production of the J/i(j state. Corresponding 
quantities for the production of ip' are shown in Fig. 4. 



3.5 r- 

3.0 

^ 2.5 - 

^2.0 - 

-^^T 1.5 - 



(a) 



J/\\f Production 



Ml Transition 
El Transition 




3.72 3.74 3.76 3.78 3.80 

Invariant mass of color-octet state (GeV) 

Fig. 3. |Vijp, -Fr, and {dn/dmf' as a function of the invariant mass mj/L's' 
of the intermediate color-octet state for the production of J/ip. 

For the production of J/ip with an Ml transition, the initial color-octet state is an L' = 
state. The results of Fig. 3 indicate that the square of the matrix element \Vij\'^ and Fr 
are approximately constants as a function of the color-octet invariant mass. The density 
of states dn/dm,j>L'S' for L' = increases with the invariant mass of the color-octet state. 
For the El transition leading to the production of J/ip , color-octet states making the El 
transition are L' = 1 states. The quantity \Vij\'^ increases monotonically from zero as the 
invariant mass increases. The factor Ff> for El transition is small because it is proportional 
to v^ (Fig. 36). Because of the centrifugal barrier, the density of states dn/dmj'L's' for L' = 1 
is zero for color-octet states in the range of invariant mass up to 3.8 GeV under consideration 
(Fig. 3c). A higher energy, L'{L' + 1)/MqR'^, is needed for the Q and Q to overcome the 
centrifugal barrier to come to a distance of i? ~ 0.5 fm for these angular momentum L' = 1 
states. As a consequence, there is no El transition from the color-octet L' = 1 state to 
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the J /if) state in the region of soft gluon emission. The Ml transition dominates over the 
El transition for J/il) production. The color-octet matrix elements for various transitions 
are given in Table II. We obtain {Oi'^(^Sf)) = O.OTeGeV^ and for {Oi^^{^S[^^)), with 
the combination {Oi^^{^S[^^)) + {3/M^){Oi^^{^S[^^)) = 0.076GeV^ which can be compared 
with the value of 0.066 extracted from the CDF measurements. 
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For the production of ijj', the size estimate gives a maximum gluon energy of the order 
of 2h/{l fm) ~ 0.4 GeV which corresponds to a maximum invariant mass of about 4.1 GeV 
for the color-octet state. For the production of ip', the radial wave function of the final state 
has a node and a change of the sign of the wave function. As a result, the matrix element 
Vij for the Ml transition involves a high degree of cancellation and the magnitude of the 
matrix element is small for the Ml transition (Fig. 4a). For the Ml transition the Fr factor 
is approximately a constant but {dn/dmjiL'S'Y increases monotonically with mj/i's' (Fig 
4). For the El transition |l^jP depends on the invariant mass of the color-octet state and 
is much larger than that for the Ml transition at mjiysi = 3.85 GeV. On the other hand, 
the Fn factor is proportional to k'^/Mg for the El transition and is small in magnitude 
(Fig 46). The square of the density of the L' = 1 states, {dn/dmy, increases with the 
color-octet invariant mass after a threshold. The combined result is a small El transition 
probability for the production oiip' as shown in Table II. We obtain (Og (^^o)) = 0.015 GeV^ 
for the Ml transition, and {Of {^ Pq)) / M^ = O.OOOlGeV^ for the El transition. The Ml 
transition probability is again much greater than the El transition probability. The sum 
(Og (^5'o)) -|- ^{Og {^Pq))/Mq agrees approximately with the matrix element extracted in 
CDF measurements H. 
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Table II. Color-Octet Matrix Elements from Eq. (^ 



Matrix Element 


Matrix Elements from 
CDF Measurements 


Matrix Elements 
Calculated with Eq. ^5$) 


{Oil^i'P,))/Ml 




0.076 GeV^ 



{Oi'^i}S,))+?>{Oi'^{'P,))/Ml 


0.066 GeV^ 


0.076 GeV=^ 


iOfeSo)) 
{Ot\'Po))lMl 




0.015 GeV^ 
0.0001 GeV^ 


{On'So))+3{Ot'('Po))/Ml 


0.018 GeV^ 


0.015 GeV^ 



It is easy to show that the probability for the transition ^Sq — > {^Sl , Sz) is the same 
for Sz = —1, 0, and 1. Therefore, an initial ^Sq state will lead to an equal population of the 
final magnetic substates of ^Si quarkonium and consequently gives an isotropic angular 
distribution of muons from the decay of the quarkonium. 

In J/ip and ip' production at fixed-target energies, the fusion of two gluons is the dominant 
mechanism of direct J/-ip and ip' production and gives rise to the color-octet ^Sq and ^PJ 
states. Thus, the direct production of J/ip and ■?/'' bound states by the emission of very soft 
gluons occurs mainly through the Ml ^5*0 — > ^5*} transition, and the produced ^5*1 bound 
state will be nearly unpolarized, with an approximately isotropic angular distribution for the 
decay muons. Assuming {OJ{^Pj))/Mq = and including the small contributions from 
the color-singlet mechanism and qq annihilation, Beneke et al. |^ obtain A = 0.15 for the 
muon angular distribution 1 + A cos^ 9 for ip' decay in the Gottfried- Jackson frame, which 
falls within the error of the experimental value of A = 0.02±0.14 for the ip' at y/s = 21.8 GeV 
T^ . The experimental A value for direct J/ip production is not yet available. The measured 



A value for total J/ip production is 0.028±0.004 at ^i = 15.3 GeV [|l|] and is -0.02±0.06 for 
y/s = 21.8 GeV |T^, which includes direct and indirect J/ip production from the decay of 
X states. The theoretical calculation of A for the total J/ip yield is still incomplete because 
of the difficulty of the Color-Octet Model to produce the observed X1/X2 ratio 0. 

We now return to the question of the universality of the matrix elements. With 
{Og^^{^So)) » {Og^^{^Pj))/ml, the discrepancy between the matrix elements in the CDF 
measurement and the fixed-target measurements is reduced to a factor of 2 for J/tp and 4 
for ip' 1^. The discrepancy can be further reduced when one takes into account the physical 
masses and allows for the finite energy carried by the soft gluon, as suggested by Beneke 
et al. p. 



VII. CONCLUSIONS AND DISCUSSIONS 

For the Color-Octet Model to be a valid description, we need to understand the soft 
gluon emission process. As the soft gluon emission occurs at a late stage of the production 
process, the final state can be described in terms of an interacting Q-Q potential. We use a 
potential model to describe the initial and final states in the soft gluon emission process. 

In our discussions of the long-distance behavior of the Q-Q potential, it is necessary to 
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take into account the effects of the spontaneous production of hght quarks as the heavy 
Q and Q pull apart. Accordingly, the potential at large distances is screened due to the 
breaking of the string. A potential which has such a property was used successfully to 
describe the color-singlet charmonium bound states and resonances. 

Not much is known about the potential between Q and Q m a, color-octet state except 
that it is repulsive at short distances. The combined potential may have a pocket deep 
enough to hold a bound state. A bound color-octet state is a relatively long-lived object 
having interesting experimental signatures such as the occurrence of chains of pions in the 
rapidity gap between the produced quarkonium state and its complementary partner as the 
octet is pulling apart from its complementary color partner. 

In the absence of evidence for the occurrence of color-octet bound states at present, we 
consider color-octet states in the continuum. For our studies, we assume a weak color-octet 
interaction which contains only the repulsive color- Yukawa interaction at short distances to 
give illustrative results concerning the production process. 

We found that the production cross section depends on three factors: the magnitude of 
the Ml and El matrix elements, the wave function or its derivatives at the origin, and the 
density of the color-octet states in the soft-gluon emitting region. It turns out that for J /ip 
production, the magnitude of the Ml matrix elements is relatively constant for the range 
of soft-gluon emissions, while the magnitude of El transition matrix elements increases 
monotonically from zero as the energy of the soft gluon increases. The El transition is, 
however, highly inhibited because within the range of low energy soft-gluons (up to uog ~ 0.8 
GeV), the density of the initial color-octet L' = 1 states is zero as a higher energy is needed 
to overcome the centrifugal barrier. As a consequence, the Ml transition dominates over 
the El transition. 

For Ip' production, the transition matrix element for the Ml transition is small because 
of the cancellation involving a final state whose wave function changes sign in different 
regions of r. However, the density of L' = 1 initial color-octet states is small and limits 
the contribution of the El transition. The Ml transition is again the dominant mode of 
transition. 

It is clear from the above discussions that the dominance of the Ml over the El transition 
arises mainly from the limiting factor of the density of P-wave color-octet states. Thus, the 
dominance of Ml soft gluon emission over the El transition is a rather general property which 
is not affected by the strength of the Q-Q potential. On the other hand, the magnitude of 
the Ml cross section is sensitive to the maximum soft gluon energy and the Q-Q interaction. 

In conclusion, while there are questions concerning the Color-Octet Model, some of 
these may be resolved by careful refinements of the details of the model. In particular, 
the dominance of the Ml transition ^Sq -^ ^S*} over the El transition ^PJ -^ ^S*} 
can explain the isotropic muon angular distribution and reduce the discrepancies of the 
matrix elements between the CDF and the fixed-target measurements. There may be the 
universality of the color-octet matrix elements when one takes into account the physical 
masses and the finite energies of the soft gluon. 
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